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f^*) ■ Abstract. 

0^ , The nonadiabatic transition probabilities in the two-level systems are calculated 

' analytically by using the monodromy matrix determining the global feature of the 

. underlying differential equation. We study the time-dependent 2x2 Hamiltonian with 

*^ the tanh-type plus sech-type energy difference and with constant off-diagonal elements 

as an example to show the efficiency of the monodromy approach. The application of 
this method to multi-level systems is also discussed. 

cr. 



1. Introduction 



Analytic calculation of the time evolution in two-level systems has been studied 
by a number of authors for a long time since the beginning years of quantum 
mechanics [Ill2lElEll3inilZllHliaiinilIIlliaiinilIll- These results have been applied to 
various areas of physics including quantum optics, laser spectroscopy, nuclear magnetic 
resonance and atomic collisions The importance of the study of quantum 

time-evolution is still increasing even now; For example, much attention has been 
paid recently to the quantum manipulation of qubits [TH] and magnetization process 
of magnetic molecules with large spin [20] . Recent rapid development of computers has 
enabled massive numerical simulation of quantum dynamics. Nevertheless, it remains 
to be important to study analytically solvable models for the following reasons: (1) in 
some ranges of physical parameters the numerical simulation becomes too difficult, and 
(2) analytic solutions give a clearer description about parameter dependence. 

Analytic solutions of quantum dynamics can be classified into several classes. Some 
of them are obtained by using hypergeometric functions. This was first found by Rosen 
and Zener jSj, which has then been generalized by several authors |7J |HJ 13 EE EH 02] • 
In these studies, the time-variable t is generally transformed into another real variable 
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z = z(t), which varies from to 1 monotonically. Then, the Schrodinger equation of 
two-level systems can be reduced to the hypergeometric differential equation, and the 
transition probability can be related to the connection problem between two pairs of 
fundamental solutions around z = and z — 1. 

One exception is the approach by Carroll and Hioe ^3]- They have studied two 
solvable classes, and in one of them they have introduced a new variable z(t) changing 
from — oo to oo as t increases and have reduced the Schrodinger equation to the Riemann- 
Papperitz equation. Recently, Ishkhanyan has pointed out that the Carroll-Hioe model 
can be understood in terms of the hypergeometric functions by considering a complex- 
valued path z(t) = (y(t) + i)/2i where y(t) is a real variable [Hj. By this complex- valued 
path, Ishkhanyan also found a new solvable class, but he did not obtain results for the 
transition probability. 

In this paper, we show that for the complex-valued path, the transition probability 
can be calculated efficiently from the 'monodromy' matrices of the corresponding 
differential equations. Monodromy is one of the global properties of differential 
equations, and has attracted much attention by mathematicians, for example, through 
the deep connection with the Painleve equations [21]. Hence, we expect that the 
monodromy approach is valuable not only because it enables one to calculate the 
transition probability for various models but also because it establishes a connection 
between physical phenomena and global features of differential equations. 

In this paper, as a concrete example, we mainly consider the following time- 
dependent two-level Schrodinger equation and obtain the transition probability using 
the monodromy associated with the solution: 



Here, the coefficients, Eq, E\ and Vo, are assumed to be real constants. This is one of 
the solvable class reported by Ishkhanyan ^3]. However, the transition probability for 
the model has not been obtained. It should be noted that this model is equivalent to 
the Rosen-Zener model [3] in the case Ex = i,, and that it also includes the special 
case of the second Demkov-Kunike model [Zj in the case Eq = 0. Hence, this model can 
give a smooth connection between the two known results. 

The plan of the paper is as follows. We give the relation between the transition 
probability and the monodromy of the hypergeometric function in §0 and the transition 
probability is calculated explicitly in § El The extension to the multi-level problems is 

| In the original Rosen-Zener model, e(t) is a constant, while V(t) has a sech-type pulse form. This 
Hamiltonian, however, is reformed by a proper unitary transformation of the wave function to coincide 
with the present model. 




(1) 



where the matrix elements are given by 



e(t) = E sech (t/T) + Ex tanh(t/T) 
V(t) = V . 



(2) 
(3) 



Analytic calculation of nonadiabatic transition probabilities 



3 



addressed in §EJ Finally, the results are summarized in §0 In Appendix A[ we describe 



the generalization of the present model and its relationship to the Carroll-Hioe's model. 
2. Hypergeometric function and monodromy 

The diagonal elements in the model ([TJ are eliminated by the following change of 
variables: 

Ci = ai exp ^i J edt^j , (4) 

c 2 = a 2 exp ^— i J edt^j . (5) 
Then, the Schrodinger equation is expressed as 

icit = V exp J q £ dtj c 2 , (6) 

ic 2t = V exp (-2i J* ed?j a- (7) 

By combining these two equations, we obtain the second-order differential equations for 
Ci and c 2 respectively as 

V t \ n 



d« + ( -2ie(t) - yj en + V 2 a = 0, (8) 



c 2ti + f 2ie(t) - c 2t + V 2 c 2 = 0. (9) 

It should be noted that the equation for c 2 is obtained by replacing e(t) by —e(i) in (jHJ). 
Hence, once the solution of the equation for a is obtained, the solution for c 2 is easily 
obtained by reversing the sign of the parameters in e(t). 

The above discussion is general. Now, we consider the specific model given by (J2J) 
and (jHJ). By substituting these specific forms of e(t) and V(t) into (jHJ) and adopting the 
change of variable as 

m = *j^2±i , (10 ) 

equation (JHJ) can be reduced to the differential equation of the hypergeometric 
function 



z(l - z)c lzz + (7 - (1 + a + P)z)c u - aPa = 0. (11) 
Here, the parameters, a, (3 and 7, are determined as 



a = iT (_ El + J E 2 + V 2^ (12 ) 



(3 = iT{-E x -sjEl + Vi), (13) 

7 = - + E T - \E\T. (14) 

In the same way, equation (0) is reduced to the hypergeometric differential equation 
with the parameters 



a' = iT(E l + JEf + V 2 ) } (15) 
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(16) 
(17) 



As already mentioned, these parameters are obtained by replacing Eq and E\ by —Eq 
and —Ei respectively in (fT^j) - (|T^ . In the following calculation, related to the calculation 
for c 2 , the prime indicates that it is obtained by reversing the sign of E and E\ from 
the original quantity without the prime. 

From (fTT)|) . it can be easily seen that the variable \z\ — > oo as \t\ — > oo. Hence, for 
discussion about the initial state it is convenient to consider the fundamental solutions 
around z = oo as 

ci = A 1 f 00 (z;a)+A 2 f 00 (z;(3), (18) 
c 2 = B 1 f O0 {z-a') + B 2 f oo {z-l3'), (19) 
where foo(z; a) and foo( z ', (3) are expressed in terms of the hypergeometric functions as 
/«,(*; a) = z- a F(a, a- 7 + l,a-/5 + l; 1/z), (20) 

(21) 



/«,(*; (3) = z^F{(3 - 7 + 1, f3, 13 - a + 1; 1/z), 



with similar definitions for foa{z;a') and / oc (2,/? / ). Since a and /3 are pure-imaginary 
in the present model, we have to choose arg(z) to determine the branch. In this paper, 
we choose 

[ — 7T/2 (£ — > +OOj 

In order to decide the initial state, it is sufficient to study the limit \z\ — > oo, in which 
case we obtain 

ci -> A lZ - a + A 2 z~ p , (23) 
c 2 -> +5 2 z" /3 '. (24) 

From (fTUj) . tpHjl and Aj's and 5;'s are determined. 

To obtain the transition probability, we assume that the initial state is the ground 
state of the Hamiltonian in the limit t — > — oo, 

-(**)• 

It may be noted that the off-diagonal elements in (|25|) do not vanish. Consequently the 
ground state wavefunction does not correspond to |oi| = 1 and a 2 = as it appears 
in the usual models. The time-evolution of the ground-state wave function is obtained 
generally as 

\ a 2 {t) ) \ -A' ) 6 



(26) 



where 



.4 



E x + JE( + V 2 



A' 



-E 1 + JE{ + V 2 



\ 2 y E i 



K 2 



(27) 
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The solution (|26|) includes an arbitary phase factor if, which is chosen zero in this paper. 
From (|2SJ), the time-evolutions of C\ and c 2 in the limit t — > — oo can be easily evaluated 

as 



E x -y/E\+V$ )t-i0o-i0l 



d(t) = ai(t)exp (i jf edtj -»• V 1 v 1 V (28) 

c 2 (t) = a 2 (t) exp (-i^edt) - -A'e" 1 ^- V^)«**>** (29) 
where the phase factors are given as 

O = T£ O |, 1 = TE 1 log2. (30) 

By using (JTUJ) and by comparing (EH"® with (j23 j) - (|2ll|l . the constants, Aj's and ^'s 
are determined as 

A a = ^ e W2-i0i-i<Ao-i^i ; A 2 = 0, (31) 

Bl = -A'e i7Ta ' /2+ ^ 1+i ^- ilp '\ 5 2 = 0, (32) 

where the phase factors, (p x and ip[, are given as 

= 2a log 2, i^ = 2a' log 2, (33) 

though these are not relevant to the calculation of the transition probability. 

By the choice of the initial condition, the time-evolution has been described only 
by the fundamental solution f^z^a) around z = oo. To be more accurate, around 
z = ioo + 1/2 (corresponding to t — > — oo) denoted by the point P in Fig. [T] (a), the 
solution is given by 

Cl (z) =A 1 f 00 {z;a), (34) 
c 2 {z)=B l f 00 {z-a l ). (35) 

On the other hand, the final state is given by the solution of the hypergeometric 
differential equation around z = —ioo + 1/2 (corresponding to t — > oo) denoted by 
the point Q in Fig. (a). The path z(t) in the complex plane is also drawn in Fig. ^ 
Then, the solution around the point Q analytically continued from the point P does 
not equal to (jMj) and (|33|) ; The solution is expressed as linear combinations of the 
fundamental solutions around z = oo. This is crucial to the calculation of the transition 
probability. 

In order to make the situation clearer, let us deform the path of z as shown in 
Fig. d (b). In this deformed path, the analytic continuation of the solution is divided 
into two parts, Ci and C 2 . Here, the path Ci denotes a round trip to the singular point 
at z = 1, while the path C 2 is a half round trip around z = oo in the clockwise direction. 
The analytic continuation along the path C 2 is easily performed, and determined only by 
the fundamental solutions around z = oo. On the other hand, the analytic continuation 
along the path Ci is nontrivial, and determined by the global character of the differential 
equation called the 'monodromy'. The monodromy is expressed by the monodromy 
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Figure 1. The path along which the analytic continuation is performed; (a) The 
original path and (b) deformed path. 



matrices as 

T(Ci) (/«,(*; «), /?)) = (/«,(«; a), /«,(*; (5)) R, (36) 
7 (C a ) (/«,(*; a), /«,(*; /?)) = (/«,(*; «'), /«,(*; /?')) (37) 

where 7(Ci) denotes the analytic continuation along the path Ci. Denoting the matrix 
elements of Q and Q' as 

*-(:;)• *-(")■ <*> 

the solutions around the point Q can be expressed by 

7(d) ( Cl (») = 1 (C 1 )(A 1 f 00 (z;a)) 

= A 1 af O0 (z;a)+Aicf oa (z;f3), (39) 
7(d) (c 2 (z)) = 1 (C 1 )(B 1 f 00 (z;a')) 

= B 1 a'f <x {z; a') + B 1 c'f 00 {z ] f3'). (40) 

Here, as shown below, the first (second) term in the final entries of ()39|) and (|40|) 
corresponds to the excited (ground) state in the limit t — > oo. Hence, in the calculation 
of the transition probability, only the element a (a') is relevant. This matrix element is 
calculated explicitly in the next section. 

Let us end this section by deriving the formula for the transition probability, using 
the monodromy matrix elements. In the limit t — > oo, the equations, ()39|) and (|40jh are 
evaluated as 

Cl{t) ^ aAie i^ a /2 e -i^v^y + . . . t (41) 

c 2 (t) - aS ie V + W/ 2e -<- El+ ^W)* + . . . (42) 

Here, we have suppressed the second term corresponding to the ground state. From 
these equations, the components of the wave function ty(t) = (ai(t), a2(t)) T is obtained 
in the limit t — > oo as 



ai(t) — ci(t)exp|— i J e(t)dt\ 
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— > aAi e i¥,+i7ra+i ^ 1_i ^° e~ 1 ^ E i +v o * + 
a 2 (t) = c 2 (t)exp j+i^ e(t)dt| 



— > a' 5 ie V+W-i^i+i0o e -i-v/^+V^* + . . . . 

By substituting (|3*T]) and (|52*)l. the wave function \l/(t) is evaluated as 

I -a'A / e 2i * 0+iW 



7 

(43) 
(44) 
(45) 



On the other hand, the Hamiltonian and the wave function of the excited state ^/e.s. i n 
the limit t — > oo are given as 



Vb -Si 



E.S. 



.4 
A' 



(46) 



where A and A' are given by (|27p. Then, the transition probability is calculated as 



P = |tf (t -> oo) ■ ^e.s. 



aA z e 



2 i-n-a— 2i0o 



/ / i/2^i7ra'+2i^ 



a A e 



2^ + t; 2 
1 

h 5 



£1 



'- -ieo _|_ a'e _ei ~V^i+^'e ieo 



oe 



V -£l ~\/ e i +1; V 



a e 



(47) 



where in the final equation we have introduced new variables, 

s = 7iTE , e x = ttTEx, v = irTV . (48) 

The last equation (|4T|) can be used for the practical evaluation of the transition 
probability. The remaining task is to calculate the elements of the monodromy matrices. 



3. Calculation of the transition probability 

One may identify several ways to calculate the monodromy matrices of the 
hypergeometric differential equations |21j . Here, we briefly explain the simplest way. 

To determine the monodromy matrix, it is crucial to use the integral representation 
of the hypergeometric function. By defining the integral 

F pq (z) = r&tt^il-ty-P-^z-t)-*, (49) 

the following relations hold: 

Pioo = c loo f (z; 0), F 0z = c 0z f (z; 1 - 7), (50) 

Fooo = Cooo/i(z;0), F lz = cuhi&l ~ ot -/?), (51) 

F i = c i/ oo (2;; a), F zoo = c 200 / 00 (z; p), (52) 

where / , fi and denote the fundamental solutions of the hypergeometric differential 
equations around z = 0,1, 00, respectively. Here, the constants, c pq 's, depend only on 



Analytic calculation of nonadiabatic transition probabilities 



8 



a, (3 and 7, and their explicit expressions are irrelevant to the present calculation. By 
applying the Cauchy's theorem to the integral in eq. ()49j). the following linear relations 
may be identified: 

F 01 + F loo + = 0, (53) 

F01 - F 0z + F lz = 0, (54) 

e(P - 7 + l)F loo - F lz - e{-a)F zoo = 0, (55) 

e{a - 7)^000 + F 0z + F ZOD = 0, (56) 
where e(-) = exp(27ri-). By eliminating F loo and F 0z , we obtain 

{F01, F zoo ) = (Fooq, Fi z )S, (57) 
S 



(5f 



e(-a) - eifl-i) 
M/3-7)-e(- 7 ) e(ar + /3-2 T )-e(/3- 7 ) \ 

V l-e(-a) e(/?- T )-l J' 

On the other hand, since the solution pair (Fooo, Fi z ) is related to the fundamental 
solutions around z = 1, the monodromy matrix for this pair is easily obtained as 

7(00(^000,^) = (Fooo, F u )r, T=(J e(7 _° a _ /?) ). (59) 

Combining (|58|l and the monodromy matrix for (.F00O5 -^Lz) is given as 

7 (C 1 )(F 01 ,F 2OO ) = (F 01 ,F ZOO )R (60) 

R=[\ l)=S->TS. (61) 

This result is easily related to the fundamental solutions, foo{z\ a) and foo(z; 0), by (|52j) 

as 

7(Ci) f/oo(z,a),— /«,(*, = f/oo(^«),— foo(z,P))R. (62) 

V coi / V coi / 

By comparing the above with the original monodromy matrix (}36|) along with ((HHJ), we 
finally obtain a = a (d = d). So, as far as a is concerned, we only need to calculate R. 
This can be performed by straightforward but slightly lengthy calculation. As a result, 
we obtain the matrix element a as 

a = e ^~7) ~e(-7) +e(-a) - 1 ^ 
e(/3 - 7) - e(a - 7) 

The matrix element a' for the solution C2 is easily obtained by reversing the sign of Eq 
and Ei in the result (|63j) . From the result for a and a', the transition probability P is 
obtained from (|4~7jl as 

„ smh^TrTiM cos 2 (ttT^ ) cosh 2 OnTEi) sin 2 (vrT£ ) , , 

P = ; 1 ; . (64) 

smh 2 {nTJEf + V 2 ) cosh 2 ! (vrT 'J Ef + V 2 ) 
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Figure 2. (a) The -Eo-dependence of the transition probability P for ttTE\ = 1. 
From the top to the bottom, ttVqT is taken 0.2,0.5,1.0,2.0. (b) The E\ -dependence 
off. From the top to the bottom, irV T is taken 0.2,0.5,1.0,2.0. The solid(dashed) 
lines show the minimum(maximum) value of P at each E\. 



Let us discuss the nature of this result. The transition probability oscillates as 
the sech-form pulse area, uTEq, changes; the transition probability has minimum and 
maximum values function of E as 

sinhVTgi) 
"mm = ; , for 7iTE = nir, 



sinh i (7rT v /E 1 2 + V c 



„ cosh (ttTEx) „ m ^ 

P max = \ U , for rrTEo = (n + 1/2)tt, 



(65) 

for 7rTE = (n+ l/2)vr, 

cosh^TrT^/^ 2 + V?) 

where n is an integer. The oscillation behavior of P for E is drawn in Fig. El (a). 

The amplitude of this oscillation becomes small as E\ increases. This feature is 
shown in Fig. |21 (b). In the case E{T 3> max(VoT, 1), we obtain the ordinary Landau- 
Zener formula 

P = e -*V *T/2 El (66) 

independent of E . 

Finally, we show that the results in the limiting cases coincide with the known 
results. In the limit Ex — > 0, the transition probability is given as 

P = g^Zgjgg) m 

cosh 2 (7rTK,)' 

which corresponds to the Rosen-Zener formula [3]. In the limit Eq — > 0, the transition 
probability is given as 

sinh 2 (7rTEi) , , 

(68) 



sinh^vrTJ^ + l/o 2 ) 



Analytic calculation of nonadiabatic transition probabilities 



10 



In this case, the present model is related to the second model in Demkov and Kunike's 
paper [7j, which corresponds to the form 

e(t) = a + 6tanh(t/T), (69) 

V(t) = c. (70) 

Their result for a = corresponds to the result (JSHJl- 

4. Application of monodromy to multi-level problems 

The application of the monodromy matrix to the transition probability is not restricted 
to the hypergeometric functions. The monodromy approach is also applicable to the 
differential equations whose monodromy is known. To show such an example, we 
consider the multi-level problem. We expect that more solvable classes can be found by 
using the present approach. 

In this section, we treat the following time-dependent Hamiltonian: 

f e(t) (i = j = l) 

V 3 (i = 1 and 2 < j < N) 

Vi (j = 1 and 2 < % < N) ' 1 ' 

(otherwise) 

where the time-dependent part e(t) is given as 

e(t) = £itanh(t/T), (72) 

and V/s (2 < j < N) are constants. It should be noted that in the limit E{T — > 
oo, this model is reduced to the extended Landau-Zener model studied by several 
authors EH 123 EH|- To eliminate the diagonal element of the Hamiltonian the 
wave function denoted by \I/(t) = (a 1; a 2 , • • ■ , dAr) T is transformed into new variables as 



a,\ exp J edt^j (i = 1) 



(73) 



Oj (2 < i < N) 

The integral in the exponent is then calculated as 

i f edt = iEiTlog(cosht/T). (74) 
Jo 

Thus, the Schrodinger equation is obtained as 



N 



, Y vAcoslit/T) 261 ^ (i = 1) 
Tci,t={ H , (75) 

Vi (cosh t/T)- 2£l d (2<i<N) 



where 



ei = LE?iT/2, ^ = -iV^T. (76) 
By changing the time variable as z = sinh(t/T), the equations are modified as 



N 



dQ_ E^(l + ^ 2 ) £l " 1/2 9 (< = 1) (77 , 

dz ~ I J=2 ' 1 ' 

Vi(l + z 2 )- £l -V2 Cl (2<i<A^) 
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We make a further change of variables as 

, 1 + z 2)- £1 -l/2 Ci ( i = 1 ) 

^ { ' ' Aj f ei + (2 < < < n) ' (7 ' 



2 + i V 2/ z + i 
where A/s are arbitrary constants satisfying 

N 

E Ai = 1. (79) 

Consequently, we finally obtain 

ddi / 1 



i 



dz V 2 

and for 2 < i < N 



N 



£i + oK + E^> ( 8 °) 

3=2 



[Z + i) ^ = Xt + V * ~ l) ^ ~ ^ ~ ^ + \)^ dr (81) 

3= ^ 



N 



This is the Okubo equation expressed by 

(zI-C)^ = Ad, (82) 
d^ 

where / is the identity operator, C is a diagonal matrix, and A is a general matrix. This 
equation has been studied by Okubo in detail j2Zj, and it is known that this form of 
equation is convenient to study the monodromy. 

Thus, it has been shown that at least one specific model of multi-level systems can 
be reduced to the differential equation whose monodromy is known. Actual calculation 
of the transition probability needs explicit treatment of the monodromy matrices, 
and remains as a future problem. The present discussion for multi-level systems is 
preliminary, and more detailed study will be needed to clarify the efficiency of the 
monodromy approach. 



5. Summary 

We have calculated the transition probability for the Hamiltonian including the tanh- 
type plus sech-type energy difference with constant off-diagonal elements. The obtained 
result gives the natural connection between the known results, the Rosen-Zener model 
and the second Demkov-Kunike model. This model also includes the Landau-Zener 
formula in the limit of the large amplitude of the tanh-type energy difference. 

In our calculation, the monodromy of the hypergeometric functions is essential. We 
have shown that the monodromy approach is also applicable to the multi-level problems. 
We expect that the use of the monodromy in the calculation of the transition probability 
does not only helps finding more solvable models but also connects global properties of 
the differential equation with the physical phenomena. Details of calculation especially 
for the multi-level problem remain as future problems. 
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Appendix A. Solvable classes 

The model considered in the main part of this paper belongs to one solvable class called 
class 1 below. It can be given as 

1 + y z at 

V(t) = {A2) 
VI + y at 

where y(t) is an 'arbitrary' monotonically increasing function satisfying y(t) — > ±oo for 
t — > ±oo. When we adopt y{t) = sinh(t/T), we obtain (j2J) and 0- For this class, the 
Schrodinger equation can be reduced to the same hypergeometric differential equation 
(fTT|) through the change of variable z(t) = (y(t) + i)/2i [Hj. Hence, all models of 
this class give the same transition probability (JMj) . In this class, however, we have to 
define the transition probability carefully. In the limit t — > — oo (y — > — oo), the matrix 
elements become 

(A,, 
(A.4, 

Hence, the wave function of the ground state in this limit has mixed components as 
treated in § El The initial state is taken as the ground state in this limiting Hamiltonian, 
and the transition probability is defined as square of modulus of the final amplitude of 
the excited states. 

The application of the monodromy is not restricted to the class 1. As discussed 
by Ishkhanyan [T3j, as long as the complex path z{t) = (y(t) + i)/2i is used, the 
calculation by the monodromy is efficient. For example, the following solvable class 
can be considered: 

™ ^ < A ' 6 > 

This class, called here the class 2, has been first studied by Carroll and Hioe |13| . 
There, the transition probability has been calculated by solving the Riemann-Papperitz 
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equation without resorting to the monodromy. By following the Ishkhanyan's discussion, 
however, our monodromy approach is also efficient for the class 2, and gives an 
alternative method. 
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